Dark Energy and the CMB 
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We find that current Cosmic Microwave Background (CMB) anisotropy data strongly constrain the 
mean spatial curvature of the Universe to be near zero, or, equivalently, the total energy density to 
be near critical — as predicted by inflation. This result is robust to editing of data sets, and variation 
of other cosmological parameters (totaling seven, including a cosmological constant). Other lines 
of argument indicate that the energy density of non-relativistic matter is much less than critical. 
Together, these results are evidence, independent of supernovae data, for dark energy in the Universe. 
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Introduction. Cosmologists have long realized that there 
is more to the Universe than meets the eye. A wide vari- 
ety of evidence points to the existence of dark matter in 
the Universe, matter which cannot be seen, but which 
can be indirectly detected by its contribution to the 
gravitational field. As observations have improved, the 
phenomenology of the "dark sector" has become richer. 
While dark matter was originally posited to explain what 
would otherwise be excessively attractive gravity, dark 
energy explains the accelerating expansion — an appar- 
ently repulsive gravitational effect. The most well-known 
argument for this additional dark component is based on 
inferences of the luminosity distances to high-z super- 
novae ^ . The anomalously large distances indicate that 
the Universe was expanding more slowly in the past than 
it is now; i.e., the expansion rate is accelerating. Acceler- 
ation only occurs if the bulk pressure is negative, and this 
could only be due to a previously undetected component. 

Here we argue for dark energy based on another grav- 
itational effect: its influence on the mean spatial curva- 
ture. This argument |^ does not rely on the supernovae 
observations and therefore avoids the systematic uncer- 
tainties in the inferred luminosity distances. It is based 
on a lower limit to the total density, and a smaller up- 
per limit on the density of non-relativistic matter. The 
lower limit comes from measurements of the anisotropy 
of the cosmic microwave background (CMB) whose sta- 
tistical properties depend on the mean spatial curvature 
which in turn depends on the mean total density. 
We find that the CMB strongly indicates that > 0.4, 
where f2 is the ratio of the total mean density to the crit- 
ical density (that for which the mean curvature would 
be zero). Upper limits to the density of non-relativistic 
matter come from a variety of sources which quite firmly 
indicate ri„j < 0.4. 

The CMB sensitivity to curvature is due to the depen- 
dence on curvature of the angular extent of objects of 
known size, at known redshifts. CMB photons that are 
penetrating our galaxy today, were emitted from a thin 
shell at a redshift of z ^ 1100 (called the "last-scattering 



surface" ) during the transition from an ionized plasma to 
a neutral medium. The "object" of known size at known 
redshift is the sound-horizon of the plasma at the epoch 
of last-scattering. Its observational signature is the loca- 
tion of a series of peaks in the angular power spectrum 
of the CMB. 

One must be careful about using current CMB data to 
determine f2 or any other cosmological parameters for 
several reasons. First, these are very difficult experi- 
ments, and the data sets they produce have low signal-to- 
noise ratios and limited frequency ranges, complicating 
the detection of systematic errors. Use of different cali- 
bration standards further increases the risk of underesti- 
mated systematic error. To counter these problems, we 
examine the robustness of our results to editing of data 
sets, and check that the distribution of model residuals 
is consistent with the stated measurement uncertainties. 

Second, the CMB angular power spectra depend on 
a number of parameters other than the curvature. To 
some degree, a change in curvature can be mimicked by 
changes in other parameters. We therefore vary six pa- 
rameters besides the curvature, placing mild prior con- 
straints on some of these so as not to explore unrealistic 
regions of the parameter space. 

Finally, existing data are insufficient to firmly establish 
the paradigm for structure formation which we have as- 
sumed: structure grew via gravitational instability from 
primordial adiabatic perturbations. Our conclusions de- 
pend on this assumption. At present, this counts as a 
possible source of systematic error. Fortunately, future 
CMB data will verify (or refute) the paradigm and will 
also allow for the determination of f2 with greatly reduced 
model dependence [Q. 

The data. Present data are already so abundant that 
it must be compressed before it can serve as the basis 
for a multi-dimensional parameter search. Fortunately, 
all data sets have been compressed to constraints on the 
angular power spectrum, Ci = 2tt J C{d)Pi{cosd)d{cos9) 
where C{9) is the correlation function. Because of the 
tremendous reduction in the size of the data sets, this 
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data compression is called "radical compression" ||]. 
Here we use the radically compressed data from 



http : //www. cita.utoronto . ca/knox/radical .html. 

In this compilation the non-Gaussianity of the power 
spectrum uncertainties has been characterized for a num- 
ber of experiments with a lot of the weight (including all 
those plotted with large symbols in Fig. 1); assuming 
Gaussianity leads to biases 0] . 
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FIG. 1. Constraints on the angular power spectrum: those 
marked with large symbols are TOCO (fiUed triangles) Q, 
CAT (open squares) §, SK (pentagons) §, 0VR05M (open 
triangle) ||| and MS AM (filled squares) 0. The model 
curves are standard COBE- normalized CDM (dotted), the 
best-fit Q = 1 model (solid), the best-fit Q = 0.4 model 
(dot-dashed), and the best-fit = 0.2 model (dashed). The 
lower panel shows residuals of the best-fit Q — 1 model. 



The Search Method. We search over a seven-dimensional 
parameter space specified by il, ftbh^, ilcdmh^, ^A^^, t, 
Us and Cio, where fli = pi/pc and i = 6, cdm, A is for 
baryons, cold dark matter and a cosmological constant 
respectively, pc = 3Hq/{8ttG) is the critical density, r 
is the optical depth to Thomson scattering, iis is the 
power-law index of the primordial matter power spec- 
trum, and Cio serves as the normalization parameter. 
The Hubble constant, Hq = 100/ikmsec~^ Mpc~^, is a 
dependent variable in this space, due to the sum rule: 
r^A + i^b + f^cdm = ^- Note that, for specificity and 
simplicity, we have chosen the dark energy to be a cos- 
mological constant; other choices (e.g., qunitessence |1I|) 
would not significantly affect our curvature constraints. 

For each value of fl we vary the 23 other parameters 
(six cosmological and 17 calibration — one for each exper- 
iment) to find the minimum value of — Xd'^ Xp- Here 
Xd is the offset log-normal form explicitly given in Eq. 
39-43 of D, which was shown to be a good approxima- 
tion to the log of the likelihood function. Information 
from non-CMB observations is included as a prior con- 
tribution, Xp- Unless otherwise stated, we assume that 
h = 0.65 ± 0.1 (a reasonable interpretation of several 
measurements |l|]) and Qbh^ = 0.019 ± 0.003 (from |l|] 



but with a 40% increase in their uncertainty). We use 
the Levenberg-Marquardt method to find the minimum 
value of for each value of il. We stop the hunt when 
the new is within 0.1 of the old value. We tested 
this method on simulated data and recovered the correct 
results. 

The likelihood of the best-fit model, C{il) is propor- 
tional to exp(— x^/2). Ideally we would marginalize over 
the non-Jl parameters rather than maximizing over them. 
However, we note that in the limit that the likelihood 
is Gaussian, these two procedures are equivalent. More 
generally, in order for marginalization to give qualita- 
tively different answers there would have to be, with de- 
creasing ri, a very rapid increase in the volume of param- 
eter space in the non-il direction with x^'s comparable to 
the minimum x^- Inspection of the Fisher matrix leads 
us to believe this is not the case. 

The Results. Our main results are shown in Fig. 2: the 
relative likelihood (oc exp(— x^/2)) of the different values 
of il. Including all the data, the best-fit (minimum x^) 
il = 1 model is 2 X 10^ times more probable than the 
best-fit il = 0.4 model, il < 0.7 is ruled out at the 95% 
confidence level. 
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FIG. 2. Relative likelihood of fi and ~x over the degrees of 
freedom (for il= V) for different collections of data sets. 

To test the robustness of this result, we edited out 
single data sets suspected of providing the most weight. 
Most of these editings produced little change. Only the 
omission of TOCO changes things substantially, and even 
then, the best-fit SI = 1 model is 150 times more probable 
than the best-fit il ~ 0.4 model. We also edited pairs of 
data sets: for no CAT and TOCO, no MSAM and CAT, 
and no MSAM and TOCO, we find = 1 to be 120, 
2.5 X 10^ and 8 times more likely than SI — 0.4. Also 
shown, as measures of goodness-of-fit, are x^ and the 
degrees of freedom. The x^ value for the "All" case is a 
bit high, but one expects even higher ones over 8% of the 
time, so there is no strong evidence for inconsistencies 
in the data. As further indication of the robustness of 
the result, one can see from the "TOCO" panel of Fig. 
2 that it persists even when all but a single data set is 
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removed. 

For the "All" case, the best-fit f2 = 1 model has il^h^ = 
0.019, h = 0.65, nA = 0.69, r = 0.17 and n = 1.12 
and is plotted in Fig. 1. There are degeneracies among 
these parameters though and none of them is strongly 
constrained on its own. For example, an equivalently 
good fit (to just the CMB data) is given by the following 
model with no tilt or reionization: fl = 1, flijh^ ~ 0.021, 
h = 0.65, flA = 0.65, r = and n = 1. 

We also covered the ilm, Ha plane, at each point find- 
ing the minimum possible with variation of the re- 
maining 5 parameters. Figure 3 shows the Ax^ = 1,4 
and 9 contours in this plane (which, for a Gaussian, cor- 
respond to the 40%, 87% and 99% confidence regions, 
respectively) . 




FIG. 3. Likelihood contours in the f2m, f^A plane. Contours 
show Ax^ = 1, 4, 9. In green is the 95% confidence region for 
fim from cluster baryon fraction determinations. The two 
solid lines encompass the allowed 2 — a region from super 
nova data ||l|. 

Discussion. Figure 3 also shows constraints on from 
clusters. Although constraints on arise from a vari- 
ety of techniques (for reviews see [^) perhaps the most 
reliable are those based on the determination of the ratio 
of baryonic matter to dark matter in clusters of galaxies 
|^-|l9| . With the assumption that the cluster ratio is the 
mean ratio (reasonable due to the large size of the clus- 
ters) 16, 1^ , |20| , and the baryonic mean density from 
nucleosynthesis, one can constrain the range of allowable 
values of flm- Since only the baryonic intracluster gas 
is detected, the upper limits on Qm from this method 
are better understood than the lower limits. Mohr et 
al. ||l^ find, from a sample of 27 X-ray clusters, that 
(including corrections for clumping and depletion of the 
gas) < (0.32 ± 0.03)/v/V0-65. Including the Hub- 
ble constant uncertainty {h = 0.65 ± 0.1) this becomes 
ri,„ < 0.32 ± 0.05. Assuming 10% of the baryons to be 
in galaxies as opposed to the g estimated by [|5|, 

we find ^Irn — 0.29 ± 0.05. Results from observations of 



the Sunyaev-Zeldovich effect in clusters are consistent, 
though less restrictive: flm = 0.31 ±0.1 Most other 
methods (those that do not rely on the cluster baryon 
fraction) generally result in formally stronger upper lim- 
its to rim- This increases our confidence in the Mohr et 
al. rim upper limit, but we do not quote these stronger 
constraints due to our concerns that they are affected by 
systematic uncertainties that are more difficult to quan- 
tify than those in the baryon fraction method. 

There have been a number of other analyses |^l| of 
CMB anisotropy data which generally obtained weaker 
constraints on . There are technical differences be- 

tween our work and previous work: we account for the 
non-Gaussianity of the likelihood function, allow for cali- 
bration uncertainties, place "sanity" priors on the Hubble 
constant and the baryon density, and vary six parameters 
in addition to the curvature. Also, much of the strength 
of our argument comes from data reported within the last 
year. 

The verdict from the CMB is now in. It does not de- 
pend on any one, or even any two, experiments. It clearly 
points towards a fiat Universe and, together with clus- 
ter data, strongly indicates the existence of dark energy. 
These conclusions are consistent with, and independent 
of, the supernovae results. The completely different sets 
of systematic uncertainties in the two arguments further 
strengthen the case. Other constraints in the Q,n,HA 
plane were recently obtained by combining cosmic 
flow data with supernovae observations. 

We have neglected several data sets, all of which, if 
included, would only strengthen our conclusions. Two 
of these are PythonV and Viper |^^. PythonV and 
Viper together trace out a peak with centroid near I — 
200, and a significant drop in power by I — 400. They 
have not been included because of the strong correlations 
in the existing reductions of the data; a new reduction 
with all correlations specified will soon be available for 
PythonV. 

Any model without a drop in power from I = 200 to 
I = 400 has difficulties agreeing with all the data. Mod- 
els fitting this description include the adiabatic models 
considered here with fl < 0.4 and also topological defect 
models, whose breadth is a consequence of the loss of the 
coherent peak structure [26| . 

We have been concentrating on implications of the 
peak location, but the height is also of interest. With 
fixed h, it is additional evidence for low flm- The 
lower Vlmh? , the later the transition from a radiation- 
dominated Universe to a matter-dominated Universe and 
the larger the early ISW effect, which contributes in the 
region of the first peak j2^]. For flat models, the best fit 
is at Vl,n — 0.4 with = 1 four times less likely. 
Conclusions. We have shown that = 1 is strongly 
favored over f7 = 0.4. This result is interesting for two 
reasons. First, f2 = 1 is a prediction of the simplest 
models of inflation. Second, together with the constraint 
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Clm < 0.4, it is evidence for dark energy. 

The CMB can say little about the nature of the dark 
energy. A cosmological constant fits the current data, 
but then so would many of the other forms of dark en- 
ergy proposed over the past few years. Generation and 
exploration of new theoretical ideas as to the nature of 
this dark energy is clearly warranted. 

Measurements of CMB anisotropy have already deliv- 
ered on their promise to provide new clues towards an 
improved understanding of cosmological structure for- 
mation and fundamental physics. We look forward to 
greater clarification of the dark energy problem, as well as 
possibly new surprises, from improved CMB anisotropy 
measurements in the near future. 
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